Small Angle Scattering by Fractal Aggregates: A Numerical Investigation of the 
Crossover Between the Fractal Regime and the Porod Regime 
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Fractal aggregates are built on a computer using off-lattice cluster-cluster aggregation models. 
The aggregates are made of spherical particles of different sizes distributed according to a Gaussian- 
like distribution characterised by a mean ao and a standard deviation a. The wave vector dependent 
scattered intensity I{q) is computed in order to study the influence of the particle polydispersity 
on the crossover between the fractal regime and the Porod regime. It is shown that, given ao, the 
location of the crossover decreases as a increases. The dependence of qc on a can be understood 
from the evolution of the shape of the center-to-center interparticle-distance distribution function. 

PACS numbers; 61.12.-q, 61.43. Bn, 61.43. Hv 



Small-angle X-ray scaitering as well as small-angle neutron scattering have been widely used to study the structure 
of disordered systemsEl^El. In particular they have allowed to demonstrate that silica aerogels are made of connected 
psj . fractal "blobs"lj'Q. In such materials, the wave vector dependent scattered intensity I{q) exhibits two crossovers related 
T-H ' to the two characteristic lengths, mean particle diameter ap and mean blob size ^. The first ongpSeparates the low-q 
{q « ^~^) saturation, called the Guinier regimeO, from the intermediate power-law behaviouroQ, called the fractal 
regime. The second one, on which we will focus in this letter, separates the fractal regime from the high-g {q >> ^) 
q^^ behaviour, called the Porod regimet . 

In some previous studies, it was assumed that the crossover between the fractal regime and the Porod regime was 
^ always exactly located at qc = I-k j a^. ThiSj-assumption was even used by some authors to quantitatively determine the 
. average diameteitlu or the gyration radiusuu of the particles. In several other papers, an average particle radius w?f 
' extracted from a fit of the scattered intensity /(g) to analytical expressions which do not include polydispersitynfl. 
^n Such kind of analysis has been done in the case of base catalysed and neutrally reacted silica aerogels where the 
particle size polydispersity impedes the determination of ao by conventional methods such as electron micrography. 
In this letter we show that such an assumption is only valid when the standard deviation a of the particle size is small 
compared to the avei|asa|ao . ■— ■ 
G \ In previous papersll3M it has been shown that cluster- cluster aggregation modelalj can satisfactorily reproduce 
i the structure of aerogels. Furthermore, if we are only interested in short range correlations, it is sufhcienJ-to consider 
^ ' a single aggregate obtained with a simplified aggregation process: the hierarchical cluster-cluster modelllj. j— , 
O I Here for simplicity we have considered the case of chemically limited cluster-cluster aggregation (CLCA)llJ using a 
O ' three dimensional off-lattice hierarchical procedure. The hierarchical scheme is an iterative method which starts with 
^ \ a collection of iVp — 2^ particles at iteration i = and ends with a unique aggregate of particles at iteration p. At 
an intermediate iteration z, one has a collection of Nc = 2^^* independent aggregates, each of them containing iV = 2* 
particles. To proceed to the next iteration, the 2^^* aggregates are grouped into pairs and is built a new aggregate 
1 with each pair according to a specific sticking rule. As soon as it is obtained, the new aggregate is randomly disoriented 
" " ' and is stored in the collection for the next iteration. The sticking rules depend on the chosen aggregation process. In 
the CLCA case the sticking rule is as follows: A particle of one cluster and a particle of another cluster are chosen 
at random as well as a random direction in space. The two clusters are disposed such that these two particles are in 
contact with their centers aligned along a random direction. Then a test of overlap is made for the other particles. If 
an overlap is found, the trial is discarded and another choice is made for both particles and random direction. Then, 
the resulting aggregate is stored for the next iteration. In this model, the resulting fractal dimcnsiontj is equal to 
about 2. 

In this work the hierarchical procedure is initialised with a collection of iV polydisperse particles whose diameters ai 
are distributed according to a truncated Gaussian distribution. We use a standard library subroutine which generates 
a set of random variables xi of zero mean and standard deviation equal to one, distributed according to the normal 
probability law: 

1 

ff(^) = 2 (1) 
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The diameters are calculated by: 

flj = 1 -I- sxi (2) 
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where s is an input parameter and where all the x^'s leading to negative a^'s have been discarded. Then we calculate 
the average oq and the standard deviation a from the usual formulae: 
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and 



ao = -> a,; 
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By varying the parameter s, one can vary the effective standard deviation a^ffl— cr/ag) which is the only relevant 
dimensionless parameter in our problem. Note that CTe// is practically equal to s for s < 0.17, where ao — 1. In figure 
1 we show typical histograms for different values of aef / • 

In the general case (where the standard deviation may be large) one can no longer calculate the scattered intensity 
I{q) as a product o£^ form factor P{q) and the structure factor S{q). One should go back to the calculation of the 
scattered amplitudal3, which is proportional to: 

i = ^ y e*« (^~'+^)d3x (4) 

where fi refers to the center of the i-th particle and x refers to a running point inside the volume of the i-th particle 
with respect to its center. The integral inside the sum, which should be performed over the volume of the i-th particle, 
can be calculated as a function of a^, assuming isotropy and homogeneity inside the sphere, leading to: 

A = Y,e''-^'Mq) (5a) 

i 

with: 

Then, assuming a random orientation of the aggregate over the direction of q, the scattered intensity I(q) = \A\'^ can 
be written as: 

m-j:^.A;^ (6) 

where = \ fi - r*, |. 

Note that, since the i and j dependent product AiAj appears inside the sum, the result cannot be split in two 
parts. In particular one cannot use the distance distribution function f(r) = ^^i) to calculate an intermediate 

structure factor S{q). Here the double sum has been calculated directly. 

The numerical results are reported in figure 2 where we have plotted the intensity I{q) as a function of the "reduced" 
wave vector qao- In figure 2a is shown the I{q) curve obtained for aeff=0-, in this curve we can observe the three 
regimes mentioned above. Figure 2b shows three I{q) curves corresponding to three aggregates of different aeff values. 
In this figure, in order to better see the crossover between the fractal and the Porod regime, we use the convenient 
I{q)q^ representation where the first maximum corresponds to the crossover. As indicated by the arrows on the left 
of figure 2b, the location of the crossover between the Guinier and the fractal regime is the same in all cases. Also, 
the slope of the fractal regime is unchanged. These two results mean that the particle size polydispersity does not 
affect the intra- aggregate long range particle correlations as well as the overall size of the aggregate. However, the 
location of the maximum corresponding to the crossover between the fractal and the Porod regime strongly depends 
on aeff - In the monodisperse case ((Te//=0), one has qcag = 27r, as expected, but, for increasing a^ff values, gcOo is 
pushed towards low values. On the other hand, in figure 2b, one observes how the oscillations of the Porod regime 
are more and more damped as the degree of polydispersity increases. 

In figure 3 we report the variation of qcUo as a function of (Jeff- There is first a quite slow decrease, but after a 
sigmoidal-like behaviour, gcQo reaches, for Ce// > 0.1, a net linear behaviour which can be approximated by: 
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1c = —(1 - l-6cre//) 



(7) 



The overall decrease of QcCIq when increasing polydispersity can be attributed to the fact that larger particles dominate 
the scattering. 

In order to understand the small change of regime observed near (Te//=0.1 in figure 3, we have calculated the 
center-to- center interparticle distribution function /(r) for aggregates of different degrees of polydispersity, which 
has been normalised as follows: 

oo 

J /(r) Wdr = (8) 



where N is the number of particles of the aggregate. 

In practice, to calculate /(r), we choose a given path Sr, and we calculate the number of interparticle distances 
lying between r and r + 5r. Then we divide the result by Airr^NSr. In figure 4 we compare /(r) curves with (7e//=0.1 
and with Ce// — 0.16 with the monodisperse curve ((Te//=0). The delta peak at r/oo = 1 as well as the disconiinuity 
at r/ao — 2, which has been attributed to short range interparticle correlations in a previous publicationtiJ, are 
progressively washed out when introducing polydispersity. Even if I{q) and /(r) are not directly related, it is worth 
noticing that the bump observed in figure 3 occurs at the value cre//=0.1 at which the two peaks merge into a single 
broad peak. Therefore the change of behaviour occurring near aeff—0-1 can be attributed to a modification of the 
short range correlations due to polydispersity. Note that such effects have nothing to do with the truncation of the 
Gaussian distribution that we have considered since this truncation becomes effective for (Jeff > 0.17. On the other 
hand, we can understand why the I{q) curve presents only one maximum at high-g for aeff > 0.1 (see figum-ih). 

We have considered other aggregation models: the diffusion limited cluster-cluster aggregation (DLCA)EZl'E3 and 
chemically limited particle-cluster aggregation (the so-called Eden modcO). The first model gives fractal aggregate 
with a fractal dimension D ~ 1.78, and the second model corresponds to a homogeneous aggregate of dimension equal 
to 3. In figure 5 we show the two I{q) resulting curves for a fixed value of (Jeff, in comparison with the CLCA case. 
In this figure, we see that the crossover between the fractal regime and Porod regime remains at the same qcUo value. 
This fact is due to the common short range correlations in all the considered cases, meaning that relation (7) is valid 
for a wide range aggregated systems. 

In summary, we have demonstrated that the crossover between the fractal regime and the Porod regime is very 
sensitive to the particle polydispersity. An empirical relation between the crossover wave vector qc, the particle average 
diameter oq and the standard deviation a was obtained. We have also shown that the results are almost independent 
on the aggregation mechanism. 

We acknowledge interesting discussions with M. Foret. One of us (A.H.) would like to acknowledge support from 
CONICIT (Venezuela). 
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FIG. 1. Histograms of the particle size Gaussian-like distribution for three different cxe// values. 

FIG. 2. (a) Log-log plot of I{q) versus qao for CTe//=0. (b) Log-log plot of I{q)q'^ versus qao for three different CTe// values. 
The white arrows indicate the crossover between the Guinier and the fractal regime. The black arrows indicate the crossover 
between the fractal and the Porod regime. All these curves result from an average over 32 simulations with N = 128 particles. 



FIG. 3. Plot of gcCio versus (Jeff- 

FIG. 4. Plots of /(r) versus r for three a^ff values. All these curves result from an average over 32 simulations with 
A'' = 128 particles. 

FIG. 5. Log-log plots of I{q)q'^ versus qao for three different models with (Te//=0.34. All these curves result from an average 
over 32 simulations with N = 128 particles. 
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Fig. 1 




Fig 2a 




Fig 2b 




I(q)q (arbitrary units) 
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